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We present a continuous time random walk model for the scale-invariant transport found in a self- 
organized critical rice pile [Christensen et al., Phys. Rev. Lett. 77, 107 (1996)]. From our analytical 
results it is shown that the dynamics of the experiment can be explained in terms of Levy flights 
for the grains and a long-tailed distribution of trapping times. Scaling relations for the exponents 
of these distributions are obtained. The predicted microscopic behavior is confirmed by means of a 
cellular automaton model. 



PACS numbers: 64.60.Lx,05.40.+j, 64.60. Ht 

Self-organized criticality (SOC), or the spontaneous 
emergence of scale invariance in nonequilibrium systems 
has attracted a great interest as an explanation of frac- 
tal behavior in nature B . Sandpiles rapidly became the 
paradigm of SOC, but it has been only very recently that 
they have been shown to be characterized by power-law 
distributions of avalanches BEl . On the other hand, un- 
derstanding the complex behavior of granular media is a 
challenge of fundamental physics per se. Between many 
other amazing properties, granular systems can behave 
simultaneously as solid or liquids and show a glassy dy- 
namics with extremely slow relaxations J4|. In addition, 
the transport properties found in granular systems dis- 
playing SOC H turns out to be very similar to the disper- 
sive transport taking place in amorphous semiconductors 
and polymers J6| . Finally, there are close connections be- 
tween sandpile models and interface depinning J7|. 

The experimental system that we want to analyze is 
a rice pile, built in the narrow gap between two vertical 
plates, over a quasi one-dimensional support of length L. 
Rice grains are added to the left side, where a vertical bar 
between the two plates forms a wall that keeps the grains 
inside the system. In contrast, the right side is open, and 
allows the exit of the grains out of the pile. Starting with 
an empty system, the slow addition of grains makes the 
pile to grow until the profile reaches the open boundary 
at the right. After this transient time, the pile arrives to a 
quasi-stationary state where the average slope fluctuates 
around a well defined angle of repose, and the influx of 
grains at the left equals on average the outflux at the 
right exit. At this point it has been shown by the Oslo 
group Q that the rice pile displays SOC if the shape of 
the grains is anisotropic enough to prevent the rolling of 
the grains down the slope, suppressing the inertial effects 
and enhancing the dissipation by means of the friction. 
This result, apart of being the first unambiguous evidence 
of SOC in granular media |g,|| remarks the fact that SOC 
is associated to strongly dissipative systems. 

Christensen et al. || have studied the transport prop- 



erties of individual grains through the Oslo rice pile. The 
transit time of tracer grains, defined as the time neces- 
sary for a tracer to escape from the pile, was measured 
experimentally. The results led to a power-law distribu- 
tion of transit times for long times preceded by a flat 



region, i.e., 



Ptr(T) 



constant for small T, 
1/T a for large T, 



(1) 



with a = 2.4 ± 0.2 ||. Moreover, when the system size 
is varied, the distribution verifies a finite size scaling, 



P tr {T,L) = L- v F(T/L V ), 



(2) 



with v = 1.5 ± 0.2 and v' = 1.4 ±0.2 §. The fact of 
having v ~ v' follows from the normalization condition. 
From here, the scaling of the mean transit time with 
system size was found to be "very anomalous": 



(T> ~ L". 



(3) 



The goal of this Letter is to study the microscopic 
properties of the transport of grains inside the rice pile, 
by means of a continuous time random walk model ||. 
Comparing our theory with the experimental findings, we 
can give a form for the distribution of trapping times and 
the distribution of flights in the real system. In addition 
we test our conclusions using a one dimensional cellu- 
lar automaton modeling the experiment p] that connects 
transport with avalanche dynamics. 

We are going to consider the motion of a single grain 
or particle through the profile as essentially one dimen- 
sional. This can be done because in the experiments the 
path of the particles takes place between two points, it 
starts at the top of the pile, next to the left wall, and 
ends at the rightmost extreme of the support. The rice 
grain remains at rest, trapped at position x during a ran- 
dom time interval t until some avalanche reaches it. At 
this point the grain performs an instantaneous jump, or 



flight, of random length I, after which it becomes trapped 
at x + I. Then the dynamics of a particle is described in 
terms of a continuous time random walk, fully specified 
by the distribution of trapping times tp(t) and by the 
distribution of flights <fi{l). To be precise ip(t)dt is the 
probability that the particle is trapped at a given posi- 
tion a time between t and t + dt, whereas <j){l)dl gives the 
probability that the particle jumps a distance between I 
and I + dl during an avalanche. Notice that we can as- 
sume I > always, since in the experiment no mechanism 
allows the grains to climb the profile, always decreasing 
to the right. This will be a great simplification in the 
calculations in comparison with models for diffusion in 
amorphous semiconductors [|9|. In addition, the length 
of the flight I will not be limited by the system size. 
Both variables t and I are taken as independent random 
processes. The assumption of statistical independence 
has succeeded in reproducing experimental results, for in- 
stance in Ref. |1(J] . Finally note that with the hypothesis 
of instantaneous jumps we are nothing else than fulfill- 
ing the usual condition for SOC, that is to have a slowly 
driven system with two separated time scales, where the 
motion of grains (or avalanches) happens at infinite ve- 
locity in comparison with any external time scale. 

The magnitude of interest is the distribution of transit 
times Ptr(T, £), where Pt r (T, C)dT gives the probability 
that the particle takes a time between T and T + dT to 
travel from the origin to position x = £ < L. In the 
context of stochastic processes the transit time is in fact 
the first passage time to level C. This distribution can 
be easily related with p(x,t), defined in such a way that 
p(x, i)dx is the probability that at time t a particle is in a 
position between x and x + dx, the time being measured 
since the addition of the particle. The probability of 
being at x < C at time t is equal to the probability of 
having a transit time T > t; mathematically, 



p(x, t)dx = 1 



P tr {T,C)dT. 



(4) 



On the other hand p(x, t) depends on the renewal density 
h(x,t), where h(x,t)dx gives the probability of jump per 
unit time in a position between x and x + dx pT| . We can 
write close equations relating p(x,t) and h(x,t), using 

ijj(t) and 4>(l): 

h(x,t) = [n<t){x) + (1 - fj,)6(x)]il>(t) 

+ [[ h{x',r)(j)(x-x')il>{t-T)dx'dr, (5) 

Jo Jo 

+ / / h(x',T)<j){x-x')^(t-T)dTdx' } (6) 

Jo Jo 

where $(i) = J. ip(t')dt' gives the probability that the 
particle survives a time larger than t trapped at any posi- 
tion, p, is the probability that at t = the particle is mov- 



ing. Equations (Q)-(||) contain the solution to our prob- 
lem, relating a measurable magnitude, P tr (T, C) with the 
magnitudes that define the microscopic dynamics, ip(t) 
and (f)(1). Applying the Laplace transform, defined as 
f(u>, s) = J °° dx J °° dt e~ UJX e~ st f(x 1 1), the equations be- 
come linear and straightforwardly solvable. In particular, 
Ptr turns out to be 



P tr (s,u>) 



(i-^(.s))[i-Mi-0M)] 

1 - $(a)$(w) 



(7) 



The Laplace transform of the mean transit time (T(C)) 
can be easily obtained from P tr (s,w) as 



<n«)> = - 



P tr (s,U>) 



i-m(i 



s =o w(l - 4>) 



-(t), (8) 



where the existence of (T(w)) depends directly on the 
existence of the first moment of ifi(t), (t) = J tip(t)dt = 
— (d^p/ds) s= Q. As in the experiment (T(£)) was found 
to be finite, we conclude that (t) exists as well. The 
same reasoning for (T 2 (£)), that was infinite in the ex- 
periment, implies that (t 2 ) is not defined and then we 
assume 

i/>(i) - B/t 2+/3 , when t -> oo, with < j3 < 1. (9) 

This means that in Laplace space [O 

t/i(s) ~ 1- (t)s + BF(-l-/3)s 1+/3 when s -> 0, (10) 

where T(-) is the Gamma function. Substituting this ex- 
pression into the equation for P tr (|7|) and inverting the 
Laplace transform for s we obtain 



Ptr(t,u) 



1-/1(1-0) B 



(11) 



for t — > ex). Comparing (||) and (|ll]) we see how the dis- 
tribution of transit times, that is a macroscopic quantity, 
is closely related with the distribution of trapping times, 
a microscopic magnitude. In fact, both are long-tailed 
distributions with the same exponent. Linking with the 
experimental result (|l|) we can conclude that the trapping 
times are power-law distributed with an exponent 



2 + /3 = a. 



(12) 



Notice that the tail of the transit time distribution only 
depends on the trapping time distribution, and not on 
the distribution of jumps. Alternatively, going back to 
Eq. (|7J) we can perform first the long distance limit, i.e., 
Cij — > 0. If we consider that the jump distribution has 
mean value but not second moment, that is, 

(f>(l) - C/l 2+1 when I -> oo, with < 7 < 1, (13) 



then <f)(u)) verifies an equation similar to (ElJ), that sub- 
stituting on (M) and inverting the Laplace transform gives 



Ptr{s,C) 



C 

1 + 7 



i-V\ 



£l+ 7 ■■ 



(14) 



(3 = 1, 



(19) 



when C — ► oo. This behavior will correspond to times 
"smaller" than C. For small times, the movement of the 
grains will be shallow and the probability of a given tran- 
sit time up to a position L will be independent on system 
size; in other words, Eq. (H4) will not depend on L. If 
we make C — L, this means that the scaling with system 
size is the same than with position. Denoting the latter 
by subindices v x andz/, i.e., P tr (T,C) = C-<T(T/C U -), 
this implies that v' = v' x , and comparing with the scaling 
found in the experiment (H) we have 



1 + 7 = 1/'. 



(15) 



One can consider other asymptotic forms for (f)(1), but 
it is only the one given by Eq. (|l|) that reproduces an 
exponent v' between 1 and 2. From here we can deduce 
that in the rice pile the scaling of Pt r given by (g) means 
that the distribution of jumps has a finite mean (I) but an 
infinite variance. This kind of distributions correspond 
to Levy flights [[L3| , and give rise to a superdiffusive be- 
havior, as one can verify by finding p(u>,s) and from here 
obtain (similarly as in Eq. (JsJ) ) that (x(t)) ~ t(l)/{t) and 
(x 2 (t))oi(l 2 )=oo. 



Moreover, if the limit s — > is performed in Eq. (14), 
one obtains the behavior for large times but "smaller" 
than £, that turns out to be independent on t: 



Ptr(t,£) 



C 



(l + 7 )(i)^+V 



when t < jp-C 



(16) 



This corresponds to the flat region observed in the transit 
time distribution before the power-law decay. In fact, 
the appearance of a plateau in P tr is an indication of the 
existence of (t). 

Now that the asymptotic form of (f)(1) is known, one 
can go back to Eqs. ( pd| ) and (g) to perform the long 
distance limit in order to obtain the scaling of P tr with 
C for large times and the scaling of the mean transit time: 



Ptr(t,C) 



B £ 

(0*2+/ 



when t > jJ-C, 



(T(C)) ~C(t)/(l) when£ 



oo. 



(17) 



(18) 



From here one obtains that v x = 1, in contrast with the 
value of v' x = 1 + 7. The reason to have v x ^ v' x is 
simple: the model does not show finite size scaling for all 
T and C. Indeed we have found scaling only for a region 
T < C(t)/(l) and T > C(t)/(l). By using numerical 
simulations we will see that this behavior is right and 
one can expect in the experiment (T(C)) ~ C as well. 
If we impose the continuity of Pt r at the crossover point 
T r ~ C we obtain 



that is, ip(t) and (f)(1) must have the same power-law tail. 
Employing also Eqs. ( |12| ) and (|l^) we get 



a = 1 



(20) 



This equation relates the exponent of the power-law tail 
of P tr with its scaling with system size, and it is well 
fulfilled by the experimental values, see Eqs. (|l|) and (j|). 
It would be difficult to design an experiment to mea- 
sure the distributions of trapping times and jumps, that 
are microscopic magnitudes, although in granular mate- 
rials microscopic and macroscopic scales are not so well 
separated as in the usual states of matter J4| . As an al- 
ternative to support our predictions we use the cellular 
automaton model introduced in Ref. || , which was found 
to reproduce the transport properties of grains quite well. 
We believe that similar results can be obtained for simi- 
lar models fl4|| . On a one-dimensional lattice, from x — 1 
to L, an integer variable h x gives the height of the pile at 
position x. Defining the local slope at x as z x = h x —h x +i 
the dynamics of the model is fully determined by the fol- 
lowing rules: if z x < z x \/x => z\ —* z\ + 1 (a grain 
is added); if z x > z x and x < L => Z x -\ — * Z x -\ + 
1, z x -* z x - 2, z x+ i — ► z x+1 + 1, and z% — ► rand(l, 2); if 
z L > z c L =*■ z L _i — > z L _i + 1, zl — »• Zl — 1, and z£ -> 
rand(l, 2); where all the sites have to be updated in par- 
allel and rand(l,2) means 1 or 2 at random, with equal 
probability. The external input of grains at X = 1 sets 
the time unit. 




t (trapping time) 
FIG. 1. Trapping time distribution ip(t) in a system of size 
L — 400. The two straight lines are power laws with expo- 
nents -.97 ± 0.05 and -2.20 ± 0.05. 

The results of Ref. § show that Eqs. (0)-(|) are 
still valid, but the exponents are determined with more 
accuracy. We reanalyze these results to obtain a — 
2.21 ± 0.05, v = 1.25 ± 0.10, and 1/ = 1.25 ± 0.10, in con- 
cordance with (|2C|). The distribution of trapping times 
can be obtained from simulations as the number of trap- 
pings of a given duration divided by the total number of 
trappings. The results are displayed in Fig. [l]. Indeed 
we obtain a power-law distribution for long times, where 
the exponent turns out to be 2 + (3 = 2.20 ± 0.05, in very 



good agreement with our prediction (|1J) if we compare 
with the independent measure of a. One can also mea- 
sure the distribution of flight lengths. The behavior for 
long distances corresponds indeed to a power law, whose 
exponent is 2 + 7 = 2.13 ± 0.05, see Fig. ||, where we re- 
strict the measure to flights starting at a fixed position. 
Observe from here that 1/, /3, and 7 are compatible with 
© and ©. 




/ (length of the flight) 
FIG. 2. Distribution of flight lengths cj>(l, L) for different 
system sizes starting from position x — 2. The length of the 
flight does not depend on L. The power-law behavior for large 
distances is characterized by an exponent —2.13 ± 0.05. 



The simulations allow one also to study the transit 
time to reach a certain position x = C smaller than the 
system size L. Keeping fixed L the exponents of the scal- 
ing of Ptr with C are v x = 1.0±0.1 and v' x = 1.2±0.1, as 
was predicted by Eqs. ( |l6| ) and ( p7| ) whereas the power- 
law exponent a = 2.18±0.05 is in good accord with (20). 



One can additionally measure the mean transit time as 
a function of C to verify that Eq. (pq) is fulfilled. 
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FIG. 3. Scaled distributions of transit times from x — 1 to 
x = jC in a system of size L — 400 for different values of C. 
The scaling is done for T>£ and T < C, resulting that the 
exponents are v x = 1.0 and v' x — 1.2, whereas the power-law 
exponent is 2. 18 ±0.05. Notice that the region with T slightly 
larger than C does not scale well. 

Moreover we can derive an additional scaling relation 
for the model. In Ref. (5| it was found that \ — v ~ F 
being \ the roughness exponent of the profile of the pile, 
that is, the fluctuations of the profile scale as L x . On 
the other hand, in Ref. J7J it was argued that D = 2 + x, 



with D the fractal dimension of the avalanches, i.e., the 
size of the avalanches scales as L D . As v ~ 1/, combining 
these relations with (E(J) one gets 



(21) 



a = D 



Taking D = 2.23 ± 0.03 @ this last result is in a fair 
agreement with the measured value of a. 

In summary, from analytical results and computer sim- 
ulations we present a coherent scenario for the transport 
in a self-organized critical granular system. The scale in- 
variance of the process is associated to long-tailed trap- 
ping time distributions and Levy flights of the grains. 
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